Using gene expression data and network topology to detect substantial pathways, clusters and switches during oxygen deprivation of Escherichia coli

<p>Abstract</p> <p>Background</p> <p>Biochemical investigations over the last decades have elucidated an increasingly complete image of the cellular metabolism. To derive a systems view for the regulation of the metabolism when cells adapt to environmental changes, whole genome gene expression profiles can be analysed. Moreover, utilising a network topology based on gene relationships may facilitate interpreting this vast amount of information, and extracting significant patterns within the networks.</p> <p>Results</p> <p>Interpreting expression levels as pixels with grey value intensities and network topology as relationships between pixels, allows for an image-like representation of cellular metabolism. While the topology of a regular image is a lattice grid, biological networks demonstrate scale-free architecture and thus advanced image processing methods such as wavelet transforms cannot directly be applied. In the study reported here, one-dimensional enzyme-enzyme pairs were tracked to reveal sub-graphs of a biological interaction network which showed significant adaptations to a changing environment. As a case study, the response of the hetero-fermentative bacterium <it>E. coli </it>to oxygen deprivation was investigated. With our novel method, we detected, as expected, an up-regulation in the pathways of hexose nutrients up-take and metabolism and formate fermentation. Furthermore, our approach revealed a down-regulation in iron processing as well as the up-regulation of the histidine biosynthesis pathway. The latter may reflect an adaptive response of <it>E. coli </it>against an increasingly acidic environment due to the excretion of acidic products during anaerobic growth in a batch culture.</p> <p>Conclusion</p> <p>Based on microarray expression profiling data of prokaryotic cells exposed to fundamental treatment changes, our novel technique proved to extract system changes for a rather broad spectrum of the biochemical network.</p

Discovering functional gene expression patterns in the metabolic network of Escherichia coli with wavelets transforms

BACKGROUND: Microarray technology produces gene expression data on a genomic scale for an endless variety of organisms and conditions. However, this vast amount of information needs to be extracted in a reasonable way and funneled into manageable and functionally meaningful patterns. Genes may be reasonably combined using knowledge about their interaction behaviour. On a proteomic level, biochemical research has elucidated an increasingly complete image of the metabolic architecture, especially for less complex organisms like the well studied bacterium Escherichia coli. RESULTS: We sought to discover central components of the metabolic network, regulated by the expression of associated genes under changing conditions. We mapped gene expression data from E. coli under aerobic and anaerobic conditions onto the enzymatic reaction nodes of its metabolic network. An adjacency matrix of the metabolites was created from this graph. A consecutive ones clustering method was used to obtain network clusters in the matrix. The wavelet method was applied on the adjacency matrices of these clusters to collect features for the classifier. With a feature extraction method the most discriminating features were selected. We yielded network sub-graphs from these top ranking features representing formate fermentation, in good agreement with the anaerobic response of hetero-fermentative bacteria. Furthermore, we found a switch in the starting point for NAD biosynthesis, and an adaptation of the l-aspartate metabolism, in accordance with its higher abundance under anaerobic conditions. CONCLUSION: We developed and tested a novel method, based on a combination of rationally chosen machine learning methods, to analyse gene expression data on the basis of interaction data, using a metabolic network of enzymes. As a case study, we applied our method to E. coli under oxygen deprived conditions and extracted physiologically relevant patterns that represent an adaptation of the cells to changing environmental conditions. In general, our concept may be transferred to network analyses on biological interaction data, when data for two comparable states of the associated nodes are made available

Divergent mechanisms underlie Smad4-mediated positive regulation of the three genes encoding the basement membrane component laminin-332 (laminin-5)

<p>Abstract</p> <p>Background</p> <p>Functional inactivation of the tumor suppressor Smad4 in colorectal and pancreatic carcinogenesis occurs coincident with the transition to invasive growth. Breaking the basement membrane (BM) barrier, a prerequisite for invasive growth, can be due to tumor induced proteolytic tissue remodeling or to reduced synthesis of BM molecules by incipient tumor cells. Laminin-332 (laminin-5), a heterotrimeric BM component composed of α3-, β3- and γ2-chains, has recently been identified as a target structure of Smad4 and represents the first example for expression control of an essential BM component by a tumor and invasion suppressor. Biochemically Smad4 is a transmitter of signals of the TGFβ superfamily of cytokines. We have reported previously, that Smad4 functions as a positive transcriptional regulator of constitutive and of TGFβ-induced transcription of all three genes encoding Laminin-332, LAMA3, LAMB3 and LAMC2.</p> <p>Methods</p> <p>Promoter-reporter constructs harboring 4 kb upstream regions, each of the three genes encoding Laminin-322 as well as deletion and mutations constructs were established. Promoter activities and TGFβ induction were assayed through transient transfections in Smad4-negative human cancer cells and their stable Smad4-positive derivatives. Functionally relevant binding sites were subsequently confirmed through chromatin immunoprecipitation.</p> <p>Results</p> <p>Herein, we report that Smad4 mediates transcriptional regulation through three different mechanisms, namely through Smad4 binding to a functional SBE site exclusively in the LAMA3 promoter, Smad4 binding to AP1 (and Sp1) sites presumably via interaction with AP1 family components and lastly a Smad4 impact on transcription of AP1 factors. Whereas Smad4 is essential for positive regulation of all three genes, the molecular mechanisms are significantly divergent between the LAMA3 promoter as compared to the LAMB3 and LAMC2 promoters.</p> <p>Conclusion</p> <p>We hypothesize that this divergence in modular regulation of the three promoters may lay the ground for uncoupled regulation of Laminin-332 in Smad4-deficient tumor cells in response to stromally expressed cytokines acting on budding tumor cells.</p

MiR-16-5p is frequently down-regulated in astrocytic gliomas and modulates glioma cell proliferation, apoptosis, and response to cytotoxic therapy

AIMS Aberrant expression of microRNAs (miRNAs) is frequent in various cancers including gliomas. We aimed to characterize the role of miR-16-5p as a candidate tumour suppressor miRNA in gliomas. METHODS Real-time PCR-based approaches were used for miRNA and mRNA expression profiling of glioma and non-neoplastic brain tissues as well as glioma cell lines. Protein levels were determined by Western blotting. In vitro analyses were performed following overexpression of miR-16-5p, trichostatin A treatment, and siRNA-mediated knock-down of HDAC3 in glioma cells. Effects of miR-16-5p on glioma cell viability, apoptosis and response to irradiation and temozolomide were assessed. RESULTS Expression of miR-16-5p was reduced relative to control brain tissue in isocitrate dehydrogenase (IDH)-mutant astrocytomas of World Health Organization (WHO) grades II, III, and IV, and a subset of IDH-wildtype glioblastomas WHO grade IV. MiR-16-5p expression was lower in IDH-mutant than in IDH-wildtype gliomas, and down-regulated in IDH-wildtype glioma lines. MiR-16-5p overexpression reduced expression of important cell cycle and apoptosis regulators in glioma cells, including CDK6, CDC25A, CCND3, CCNE1, WEE1, CHEK1, BCL2, and MCL1. In line, CDK6, WEE1, CHEK1, BCL2, and MCL1 transcript levels were increased in WHO grade III or IV gliomas. Trichostatin A treatment and HDAC3 knockdown in glioma cells induced miR-16-5p up-regulation and reduced expression of its targets. Moreover, miR-16-5p overexpression inhibited proliferation and induced apoptosis in various glioma cell lines and increased sensitivity of A172 glioma cells to irradiation and temozolomide. CONCLUSION Reduced expression of miR-16-5p contributes to glioma cell proliferation, survival, and resistance to cytotoxic therapy. This article is protected by copyright. All rights reserved

The landscape of viral associations in human cancers

Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, for which whole-genome and—for a subset—whole-transcriptome sequencing data from 2,658 cancers across 38 tumor types was aggregated, we systematically investigated potential viral pathogens using a consensus approach that integrated three independent pipelines. Viruses were detected in 382 genome and 68 transcriptome datasets. We found a high prevalence of known tumor-associated viruses such as Epstein–Barr virus (EBV), hepatitis B virus (HBV) and human papilloma virus (HPV; for example, HPV16 or HPV18). The study revealed significant exclusivity of HPV and driver mutations in head-and-neck cancer and the association of HPV with APOBEC mutational signatures, which suggests that impaired antiviral defense is a driving force in cervical, bladder and head-and-neck carcinoma. For HBV, HPV16, HPV18 and adeno-associated virus-2 (AAV2), viral integration was associated with local variations in genomic copy numbers. Integrations at the TERT promoter were associated with high telomerase expression evidently activating this tumor-driving process. High levels of endogenous retrovirus (ERV1) expression were linked to a worse survival outcome in patients with kidney cancer

Sumoylation of Smc5 Promotes Error-free Bypass at Damaged Replication Forks

Replication of a damaged DNA template can threaten the integrity of the genome, requiring the use of various mechanisms to tolerate DNA lesions. The Smc5/6 complex, together with the Nse2/Mms21 SUMO ligase, plays essential roles in genome stability through undefined tasks at damaged replication forks. Various subunits within the Smc5/6 complex are substrates of Nse2, but we currently do not know the role of these modifications. Here we show that sumoylation of Smc5 is targeted to its coiled-coil domain, is upregulated by replication fork damage, and participates in bypass of DNA lesions. smc5-KR mutant cells display defects in formation of sister chromatid junctions and higher translesion synthesis. Also, we provide evidence indicating that Smc5 sumoylation modulates Mph1-dependent fork regression, acting synergistically with other pathways to promote chromosome disjunction. We propose that sumoylation of Smc5 enhances physical remodeling of damaged forks, avoiding the use of a more mutagenic tolerance pathway.Ministerio de Ciencia, Innovacion y Universidades (BFU2015-71308-P, PGC2018-097796-B-I00)AGAUR-Generalitat de Catalunya (2017-SGR-569

qPCR in gastrointestinal stromal tumors: Evaluation of reference genes and expression analysis of KIT and the alternative receptor tyrosine kinases FLT3, CSF1-R, PDGFRB, MET and AXL

<p>Abstract</p> <p>Background</p> <p>Gastrointestinal stromal tumors (GIST) represent the most common mesenchymal tumors of the gastrointestinal tract. About 85% carry an activating mutation in the <it>KIT </it>or <it>PDGFRA </it>gene. Approximately 10% of GIST are so-called wild type GIST (wt-GIST) without mutations in the hot spots. In the present study we evaluated appropriate reference genes for the expression analysis of formalin-fixed, paraffin-embedded and fresh frozen samples from gastrointestinal stromal tumors. We evaluated the gene expression of <it>KIT </it>as well as of the alternative receptor tyrosine kinase genes <it>FLT3</it>, <it>CSF1-R</it>, <it>PDGFRB</it>, <it>AXL </it>and <it>MET </it>by qPCR. wt-GIST were compared to samples with mutations in <it>KIT </it>exon 9 and 11 and <it>PDGFRA </it>exon 18 in order to evaluate whether overexpression of these alternative RTK might contribute to the pathogenesis of wt-GIST.</p> <p>Results</p> <p>Gene expression variability of the pooled cDNA samples is much lower than the single reverse transcription cDNA synthesis. By combining the lowest variability values of fixed and fresh tissue, the genes <it>POLR2A</it>, <it>PPIA</it>, <it>RPLPO </it>and <it>TFRC </it>were chosen for further analysis of the GIST samples. Overexpression of <it>KIT </it>compared to the corresponding normal tissue was detected in each GIST subgroup except in GIST with <it>PDGFRA </it>exon 18 mutation. Comparing our sample groups, no significant differences in the gene expression levels of <it>FLT3</it>, <it>CSF1R </it>and <it>AXL </it>were determined. An exception was the sample group with <it>KIT </it>exon 9 mutation. A significantly reduced expression of <it>CSF1R</it>, <it>FLT3 </it>and <it>PDGFRB </it>compared to the normal tissue was detected. GIST with mutations in <it>KIT </it>exon 9 and 11 and in <it>PDGFRA </it>exon 18 showed a significant <it>PDGFRB </it>downregulation.</p> <p>Conclusions</p> <p>As the variability of expression levels for the reference genes is very high comparing fresh frozen and formalin-fixed tissue there is a strong need for validation in each tissue type. None of the alternative receptor tyrosine kinases analyzed is associated with the pathogenesis of wild-type or mutated GIST. It remains to be clarified whether an autocrine or paracrine mechanism by overexpression of receptor tyrosine kinase ligands is responsible for the tumorigenesis of wt-GIST.</p

In silico SNP analysis of the breast cancer antigen NY-BR-1

Background: Breast cancer is one of the most common malignancies with increasing incidences every year and a leading cause of death among women. Although early stage breast cancer can be effectively treated, there are limited numbers of treatment options available for patients with advanced and metastatic disease. The novel breast cancer associated antigen NY-BR-1 was identified by SEREX analysis and is expressed in the majority (>70%) of breast tumors as well as metastases, in normal breast tissue, in testis and occasionally in prostate tissue. The biological function and regulation of NY-BR-1 is up to date unknown. Methods: We performed an in silico analysis on the genetic variations of the NY-BR-1 gene using data available in public SNP databases and the tools SIFT, Polyphen and Provean to find possible functional SNPs. Additionally, we considered the allele frequency of the found damaging SNPs and also analyzed data from an in-house sequencing project of 55 breast cancer samples for recurring SNPs, recorded in dbSNP. Results: Over 2800 SNPs are recorded in the dbSNP and NHLBI ESP databases for the NY-BR-1 gene. Of these, 65 (2.07%) are synonymous SNPs, 191 (6.09%) are non-synoymous SNPs, and 2430 (77.48%) are noncoding intronic SNPs. As a result, 69 non-synoymous SNPs were predicted to be damaging by at least two, and 16 SNPs were predicted as damaging by all three of the used tools. The SNPs rs200639888, rs367841401 and rs377750885 were categorized as highly damaging by all three tools. Eight damaging SNPs are located in the ankyrin repeat domain (ANK), a domain known for its frequent involvement in protein-protein interactions. No distinctive features could be observed in the allele frequency of the analyzed SNPs. Conclusion: Considering these results we expect to gain more insights into the variations of the NY-BR-1 gene and their possible impact on giving rise to splice variants and therefore influence the function of NY-BR-1 in healthy tissue as well as in breast cancer

The landscape of viral associations in human cancers

Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, for which whole-genome and-for a subset-whole-transcriptome sequencing data from 2,658 cancers across 38 tumor types was aggregated, we systematically investigated potential viral pathogens using a consensus approach that integrated three independent pipelines. Viruses were detected in 382 genome and 68 transcriptome datasets. We found a high prevalence of known tumor-associated viruses such as Epstein-Barr virus (EBV), hepatitis B virus (HBV) and human papilloma virus (HPV; for example, HPV16 or HPV18). The study revealed significant exclusivity of HPV and driver mutations in head-and-neck cancer and the association of HPV with APOBEC mutational signatures, which suggests that impaired antiviral defense is a driving force in cervical, bladder and head-and-neck carcinoma. For HBV, HPV16, HPV18 and adeno-associated virus-2 (AAV2), viral integration was associated with local variations in genomic copy numbers. Integrations at the TERT promoter were associated with high telomerase expression evidently activating this tumor-driving process. High levels of endogenous retrovirus (ERV1) expression were linked to a worse survival outcome in patients with kidney cancer